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Abstract 

We discuss how correlations among momenta of the decay products may yield information 
about the CP-parity of the Higgs particle. These are correlations of decay planes defined 
by the momenta of pairs of particles as well as correlations between energy differences. 
Our study includes finite-width effects. In the narrow-width approximation the angular 
correlations coincide with previously reported results. The correlations between energy 
differences turn out to be a much better probe for CP determination than the previously 
suggested angular correlations, especially for massive Higgs bosons. 
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When some candidate for the Higgs particle is discovered, it becomes imperative to 
establish its properties, other than the mass. While the standard model Higgs boson is even 
under CP, alternative and more general models contain Higgs bosons for which this is not 
the case. In specific models, production cross sections and branching ratios will be different 
[1] but it would clearly be valuable if one could establish its CP on more general grounds, 
independent of the production mechanism. We shall here discuss how decay distributions 
may shed light on its intrinsic parity, in particular, how one can determine whether it is 
even or odd under CP. 

This situation is similar to the classical one of determining the parity of the n° from 
the angular correlation of the planes spanned by the momenta of the two Dalitz pairs [2, 3], 

tt° ^77^ (e+e-)(e + e-). (1) 

In that case, as pointed out by Yang [2], the correlation of the planes is unam- 
biguously determined by the parity-conserving 7r°77 and 7e + e~ couplings. Classically, 

— » — » — > — » 

the 7r u 77 coupling is proportional to E ■ B, where E and B are the electric and mag- 
netic field strengths, respectively. In relativistic notation, this corresponds to the coupling 
ie^ vpa F w F V(y , where F w is the electromagnetic field strength tensor. If we describe the 
photons by the vector field A^ 1 , then, with pi and P2 the momenta of the two photons, this 
7r°77 coupling has the form [3] 

iQesiplpDe^p^a, (2) 

where g e s is some effective coupling that depends on the kinematics. For comparison, the 
coupling of the standard-model Higgs to vector bosons (of mass my) is given by 

i (2 2^)^ m^T, (3) 

with Gp the Fermi constant. 

In the electroweak interactions, parity is known to be violated. However, we shall 
here assume that parity violation is limited to the couplings between vector bosons (W 
and Z) and fermions (quarks and leptons), as is the case in the standard model. Then 
it turns out that correlations similar to those of Dalitz pairs in 7r° decay determine the 
CP-parity of the Higgs boson. 

2 



In non-standard or extended models of the electroweak interactions, there exist 
"Higgs-like" particles having negative CP. An example of such a theory is the minimal 
super symmetric model (MSSM) [4], where there is a neutral CP-odd Higgs boson, often 
denoted A and sometimes referred to as a pseudoscalar. 

We shall consider the decay where a standard-model Higgs (H) or a 'pseudoscalar' 
Higgs particle (A) decays via two bosons (H^ + H^~ or ZZ), to two non-identical fermion- 
antifermion pairs, 



The two vector bosons need not be on mass shell. 

If we let the momenta (qi, q 2 , (fe, and g 4 ) of the two fermion-antifermion pairs (in 
the Higgs rest frame) define two planes, and denote by the angle between those two 
planes, then we shall discuss the angular distribution of the decay rate T, 

-- (5) 

r d0 K J 

and a related quantity, to be defined below, both in the case of CP-even and CP-odd Higgs 
bosons. 

Related studies have been reported by [5, 6] in the context of correlations between 
decay planes involving scalar and (technicolor) pseudoscalar Higgs bosons. It should be 
noted that the present discussion is more general, since finite-width effects of the vector 
bosons are taken into account. This enables us to investigate the angular correlations for 
Higgs masses below the threshold for decay into real vector bosons. In addition, we discuss 
correlations between energy differences. It turns out that under suitable experimental 
conditions these correlations provide an even better signal for CP determination. 

We let gy and qa denote the vector and axial-vector parts of the couplings, such 
that the fermion-fermion-vector coupling is given by 



In order to parametrize the vector and axial couplings, we define the couplings involving 
vector bosons V\ and V 2 in terms of angles Xi an d X2 as 



{H,A)^V l V 2 ^{hh){hh). 



(4) 




(6) 



gy = g% cos Xi, 



9 a = QiSinxi, 



i = 1,2. 



(7) 
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The only reference to these angles is through sin 2%. Relevant values are given in table 1 
[7]. Furthermore, g\ and g 2 contain CKM matrix elements and electric and weak isospin 
charges relevant to the fermions in question. The couplings of A and H to the vector bosons 
are given by (2) and (3), respectively. In the massless fermion approximation, we find 



d 8 r, ; = a 



Xi + sin(2xi)sin(2x 2 )^ dLips(m 2 ;gi,g 2 , 93,^4), i = H,A, (8) 



with dLips(m 2 ; q±, q 2 , 93, Q4) denoting the Lorentz-invariant phase space, m the Higgs mass, 
and with the momentum correlations given by 

= (gi -93X92 -94) + (gi -94X92 -93), (9) 

Yu = (91-93X52 -54) -(51 -54X52 -53), ( 10 ) 
X A = -2[(q l -q 2 )(q 3 -q A )} 2 -2[(q 1 -q 3 )(q 2 -q A ) -(q 1 -q 4 )(q 2 -q 3 )} 2 

+ (9i • 92X93 • 9 4 ){[(9i • 9s) + (92 • 9 4 )] 2 + [(91 • 94) + (92 • 9s)] 2 }, (11) 

Y A = (91 • 92) (93 • 94) [(91 -92) • (93 + 94)] [(93 - 94) • (91 + 92)]- (12) 

The normalizations are given as 

V2G F m^ A n #cff( s i> s 2) n , v , 1Q x 

C K = -D(s 1 ,s 2 ) and C A = D(s 1 ,s 2 ), (13) 

m 4m 



with 



and 



*>>^-Ri*,-J*+>*n N <' (14) 



Si = (9i + 92) 2 , s 2 = (q 3 + q A ) 2 . (15) 

Here, A^- is a colour factor, which is three for quarks, and one for leptons. Finally, m v and 
Ty denote the mass and total width of the relevant vector boson, respectively. 

We first turn to a discussion of angular correlations. The relative orientation of the 
two planes is defined by the angle 0, 

, (9! x 92) • (93 x q 4 ) ( s 

cos<p = — — — — — — . (16) 

9i x 92 93 x 94 
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We find 

cWh = SV2G F< D( Sl , , 2 ) VA (,^ r + i2sis2 

d0 dsids 2 9 (87r) b ?7r v 7 

-|- J y / sis 2 (m 2 - si - s 2 ) sin(2xi) sin(2x 2 ) cos + 2siS 2 cos 20 

d3FA - (gi, gg) g(gi, gg) [A(to 2 , 8l , S2 )] 3/2 _ 
d0d Sl ds 2 _ 9 (Stt)^ Sl ' 2 [4 008 ^ 

where A (x, y, z) = x 2 + y 2 + z 2 — 2 (xy + xz + yz) is the usual two-body phase space 

function. The term Y\ of eq. (12) does not contribute in eq. (18). 

In order to obtain a more inclusive distribution, we shall next integrate over the 

invariant masses of the two pairs. Thus, 

dr _ r m ' 2 r(™-V^f dg d 3 r 
d0 Jo 1 Jo 2 d<fi dsi ds 2 

Integrating over < < 2ti in the standard-model case, we confirm the result quoted in 

ref. [8] for Higgs decay into four fermions including finite-width effects. 

We introduce the ratios 



„=(^) 2 and (20) 



and change variables according to x\ = s\/m 2 and x 2 = s 2 /m 2 . This enables us to define 
the integrals 



F(m) 



dxi 



o (xi - pf + /i7 



Al-^l) 2 J\ (l,Xi,X 2 ) 

x dx 2 - x -J [A(l,xi,x 2 ) + 12xix 2 ] , (21) 

Jo (x 2 -/i) +/i7 

A( m) , /' f 1 -^ d, 2 vg^P (1 _ xi _ X2) , (22) 

Jo (xi — //J + /ry Jo (x 2 — //J + pq 

Z 1 dxixi /(W^) 2 x 2V /A(l,X!,x 2 ) 

5(to) = / -2- / dx 2 — * -2- . (23) 

Jo (xi — //J + pq Jo [x 2 - n) + pq 

The distributions of eq. (5) then take the form 

27T dr 

— — — = 1 + a (to) sin(2xi) sin(2x 2 ) cos0 + f3 (to) cos 20, (24) 

I n do 



where 

, , / 3tt\ 2 Aim) n . . B(m) . . 

am = — -=7—7, /5 m =2— t-4. 26 

v ; V 4 / F(m)' PV 7 F(m) v ; 

The functions a and f3 depend on the ratios of the masses of the vector bosons to the 
Higgs mass. They are given in fig.l, for values of m up to 1000 GeV. In the narrow- width 
approximation these decay correlations are identical to the ones reported in [5] and our 
analysis fully justifies this approximation. 

However, our analysis is valid also below the threshold for producing real vector 
bosons, m < 2 my. A representative set of angular distributions are given in fig. 2 for the 
cases H -> {l + n){bc) and H -> (66) for m = 70, 150, 300 GeV and m = 70, 300 GeV, 

respectively. (Of course, jet identification is required for this kind of analysis.) With 
being defined as the angle between two oriented planes, it can take on values < < it. 
There is seen to be a clear difference between the CP-even and the CP-odd cases. We 
observe that the distribution (24) becomes uncorrelated in the heavy Higgs limit, whereas 
the distribution (25) is independent of the Higgs mass. 

Let us now turn to a discussion of energy differences. In order to project out the 
lA-term of eq. (12), we multiply (8) by the energy differences (u>i — cj 2 )(^3 — 0Ji) before 
integrating over energies. In analogy with eq. (8), we introduce 



d 8 r, ; = a 



with 



Xi + sin(2xi) sin(2x 2 )Vi] dLips(m 2 ; q 1 , q 2 , q 3 , g 4 ), (27) 

Xi = Xi(ui - uj 2 )(uj 3 - (28) 

Yi = Y^i-waXwa-c^), (29) 

for i = H, A. The distribution takes the form 

2tt df n (rn) 

— — r— = 1 H : — 77; — r — r~rz — r COS 0, 30 

T d0 sin(2xi)sm(2x 2 ) V ' 

in the CP-even case, whereas there is no correlation in the CP-odd case; i.e. k (m) = 

and in this case the term Xj^ of eq. (28) does not contribute. Here, 

, x 1 /3tt\ 2 K(m) . . 

K(m) = I (le) 77m) ' (31) 
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with 



j, x f 1 dxi x 1 /^-v^) 2 x 2 [A(l,rri,a;2)] 3/2 f . 

J(m) = / / dx 2 — -2- , (32) 

Jo (xi - fx) + fij Jo (x 2 -ii) + //7 



^ \ _ f 1 dxt Jxj 
K[m) = / \ 

JO (Xi - jJii) + jl! 



7i 



x / dx 2 V } 1 1 \ 2 U 2)1 (l-x 1 - x 2 ) . (33) 

Jo {x 2 -fi) +/17 

The function k is given in fig. 3, for values of m up to 1000 GeV. In the narrow-width 
approximation 

, s 1 aim) 

K < m) = is m- (34) 

The correlation is significant for any Higgs mass, and in particular K(m) oc m 2 in the 
heavy-Higgs limit. A representative set of distributions (30) are given in fig.4 in the cases 
H -> (l + vi)(bc) and H -> (l + l-)(bb) for m = 70,300,500 GeV and m = 70,300 GeV, 
respectively. We see that the energy-weighted distribution is a much more sensitive probe 
for CP determination than the purely angular distribution of eq. (5). In this latter case, we 
compare an uncorrelated distribution with a strongly correlated one. Moreover, for V = Z 
the sin 2x-factors provide an enhancement in the energy- weighted correlations. 
This research has been supported by the Research Council of Norway. 



References 

[1] J.F. Gunion, H.E. Haber, G. Kane and S. Dawson, The Higgs Hunter's Guide, 
Addison- Wesley, 1990. 

[2] C.N. Yang, Phys. Rev. 77 (1950) 242, 722; 

N. Kroll and W. Wada, Phys. Rev. 98 (1955) 1355. 

[3] R.H. Dalitz, Proc. Phys. Soc. (London) A64 (1951) 667. 

[4] P. Fayet, Phys. Lett. 69B (1977) 489; 

S. Dimopoulos and H. Georgi, Nucl. Phys. B193 (1981) 150; 



7 



H.-P. Nilles, Phys. Rep. C110 (1984) 1; 

H. E. Haber and G. L. Kane, Phys. Rep. C117 (1985) 75. 

[5] C. A. Nelson, Phys. Rev. D37 (1988) 1220. 

[6] M. J. Duncan, G. L. Kane and W. W. Repko, Phys. Rev. Lett. 55 (1985) 579; 
M. J. Duncan, G. L. Kane and W. W. Repko, Nucl. Phys. B272 (1986) 517. 

[7] Particle Data Group, Review of Particle Properties, Phys. Rev. D45, 
Part 2 (June 1992). 

[8] B. A. Kniehl, Phys. Lett. B244 (1990) 537. 



8 



V 


/ 


sin 1\ 


w 


any 


1 


z 


e, /I, r 


0.1393 


z 


u, c, t 


0.6641 


z 


d, s, b 


0.9349 



Table 1: The factors sin 2% that give the relative importance of the axial couplings. 
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Figure captions 

Fig. I. The coefficients a and (3 of the angular correlations in eq. (24), for a CP-even 
Higgs of mass m. 

Fig. 2. Characteristic angular distributions of the planes defined by two Dalitz pairs for 
CP-even Higgs particles decaying via two W's and two Z's, compared with the 
corresponding distribution for a CP-odd Higgs particle (denoted A). Different 
Higgs masses are considered in the CP-even case. 

Fig. 3. The coefficient k of the energy-weighted angular correlation of eq. (30), for a 
CP-even Higgs of mass m. For a CP-odd Higgs, k — 0. 

Fig. 4. Characteristic energy-weighted distributions for CP-even Higgs particles decaying 
via two W's and two Z's, compared with the corresponding distribution for a CP- 
odd Higgs particle (A). 
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Figure 1 



11 



Figure 2 
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Figure 3 
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Figure 4 
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